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ABSTRACT 

A  high  fineness  ratio  body  of  revolution  with  an  ogive 
forebody  shape  was  tested  at  varying  angles  of  attack, 
Reynolds  numbers,  and  spin  rates  to  measure  the  Magnus 
forces  developed.   The  experimental  data  were  compared 
with  previously  derived  theoretical  equations  for  an  open- 
ended  cylinder  of  revolution  to  substantiate  the  theoretical 
work. 
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NOMENCLATURE 

CY  -  side  force  coefficient  based  on  frontal  area 

D  -  diameter  of  missile 

L  -  length  of  missile 

p  -  dimensionless  spin  rate 

q  -  freestream  dynamic  pressure 

R.  -  Reynolds  number  based  on  missile  length 

Li 

Rp  -  cross  flow  Reynolds  number  based  on  missile  diameter 

S  -  missile  maximum  cross  sectional  area 

T  -  air  temperature 

U  -  freestream  velocity 

Y  -  Magnus  side  force 

Ct  -  angle  of  attack 

r  -  circulation 

^L  -  coefficient  of  viscosity 

n  -  air  density 

CJ  -  spin  angular  velocity 


INTRODUCTION 

The  study  of  flight  trajectories  of  spin  stabilized 
missiles  requires  detailed  analyses  of  the  aerodynamic 
forces  developed  during  flight.   The  projectile's  lift, 
drag,  and  side  force  must  be  accurately  predicted  for 
studies  of  this  kind.   This  paper  investigates  the  flow 
over  a  high  fineness . ratio  spinning  body  of  revolution  to 
determine  the  side  force  and  its  center  of  pressure  gener- 
ated during  flight.   The  side  force  causes  drift  from  the 
plane  of  firing  and  so  is  an  important  aerodynamic  force. 

Determining  the  flow  over  a  missile  is  a  very  compli- 
cated three-dimensional  problem  in  which  the  primary 
mechanism  for  production  of  the  side  force  is  shed 
vorticity  caused  by  the  eventual  separation  of  the  boundary 
layer.   Because  of  this  the  flow  problem  is  dominated  by 
viscous  effects.   Several  theories  have  been  developed  to 
date  to  predict  the  side  force  and  moment  on  open-ended 
cylinders,  assuming  laminar  boundary  layers.   These 
theories,  however,  have  yet  to  be  substantiated  experi- 
mentally in  their  appropriate  regions  of  validity.   For 
this  reason  a  suitable  testing  apparatus  and  test  program 
were  devised  to  substantiate  these  theories.   This  thesis 
presents  work  accomplished  to  date  on  such  studies. 


HISTORICAL  BACKGROUND 

The  study  of  spinning  shapes  and  the  side  forces 
developed  by  them  has  been  carried  on  for  many  years.   To 
properly  understand  the  complicated  flow  over  a  high  fine- 
ness ratio  body  of  revolution  it  is  useful  to  study  past 
investigations  made  on  less  complex  spinning  bodies. 

In  1671  G.  T.  Walker   [1]   recorded  the  first  descrip- 
tion of  drift  of  a  spinning  body  in  flight,  the  curve  of 
a  sliced  tennis  ball.   Nearly  200  years  later  G.  Magnus   [2] 
published  the  results  of  his  experiments  to  account  for  the 
drift  of  a  spinning  musket  shot.   Magnus'  method  was  to 
fire  a  musket  ball  whose  center  of  mass  was  not  at  its 
geometric  center  and  thus  impose  a  spin  on  the  ball  as  it 
left  the  gun  barrel.   Magnus  found  a  clockwise  rotation 
would  cause  the  ball  to  drift  right  and  a  counter-clockwise 
rotation  would  cause  a  drift  to  the  left.   This  substanti- 
ated his  theory  that  an  aerodynamic  force  was  responsible 
for  the  drift.   Any  such  side  force  produced  by  a  spinning 
body  is  today  called  a  Magnus  force. 

Lord  Rayleigh  published  a  paper  in  1877  entitled  "On 
the  Irregular  Flight  of  a  Tennis  Ball"   [3]   in  which  he 
mathematically  described  the  Magnus  effects  in  terms  of  a 
potential  flow  around  a  cylinder  with  circulation.   Missing 
from  his  theory,  however,  were  skin  friction  effects 


between  the  fluid  and  cylinder  which  related  the  cylinder 
spinning  rate  to  the  magnitude  of  the  assumed  circulation. 

Little  or  no  theory  followed  Rayleigh's  work  until 
after  World  War  I.   Shortly  after  the  war  a  notable  attempt 
was  made  to  use  the  Magnus  effect  commercially  by  a  German 
named  Anton  Flettner   [4] .   He  replaced  the  sails  of  a  ship 
with  a  large  powered  rotor.   In  a  cross-wind  the  Magnus 
effect  would  power  the  ship.   Two  Flettner  Rotor  Ships  were 
built  and  one  sailed  across  the  Atlantic,  but  the  speed  and 
reliability  of  screw  propulsion  proved  to  be  the  rotor 
ship's  downfall. 

In  1924  Prandtl   [5]   published  his  work  on  the 
problem  of  creation  of  circulation  about  a  spinning  body. 
He  was  able  to  show  that  circulation  was  developed  and 
could  be  accounted  for  by  shed  boundary  layer  vorticity. 
The  net  circulation  was  equal  but  opposite  in  sign  to  the 
shed  vorticity. 

The  Second  World  War  provided  graphic  evidence  of 
the  need  to  accurately  predict  weapon  trajectories.   To 
place  a  shell  fired  from  a  gun  on  target  usually  required 
more  than  one  shot.   Corrections  were  made  from  previous 
impacts  to  achieve  the  desired  results.   Significant 
reduction  in  the  number  of  firings,  time,  and  cost  could 
be  made  if  a  weapon's  trajectory  could  be  accurately 
predicted. 

In  1953  Martin   [6]   published  the  first  analytical 
method  for  predicting  Magnus  forces  developed  by  spinning 


open-ended  cylinders  of  revolution  at  small  angles  of 
attack  in  incompressible  laminar  flow.   Kelly   [7]   followed 
in  1954  with  a  more  detailed  discussion  of  Martin's  theory. 
Power   [8],  in  1971,  carried  Kelly's  theory  to  a  high  order 
of  approximation  and  included  inviscid  circulation  develop- 
ment to  obtain  an  equation  for  Magnus  force  coefficient. 
He  improved  previous  theories  by  estimating  shed  vorticity 
that  lies  close  to  the  body  using  a  cross-flow  analogy. 
This  theory  best  models  the  essential  Magnus  force  contri- 
butions and  was  used  to  correlate  the  wind  tunnel  data. 

Martin,  Kelly,  and  Power  all  developed  their  theories 
by  studying  open-ended  cylinders.   Kelly  noted  in  his  work 
that  there  were  no  known  measurements  with  which  he  could 
compare  his  results  and  that  his  theory  was  not  directly 
applicable  to  bodies  with  different  nose  shapes.   Testing 
open-ended  cylinders  in  a  flow  is  impractical  and  nearly 
impossible  because  there  is  no  method  to  suspend  the  body 
in  the  flow  and  spin  it  without  complicated  f low-disturbin 
equipment.   For  this  reason  all  the  testing  that  is 
described  in  this  paper  was  done  using  a  cylinder  with  an 
ogive  nose  shape.   The  ogive  forebody  shape  was  selected 
because  of  its  constant  surface  curvature  which  allows 
theoretical  boundary  layer  calculations  to  be  made  more 
easily. 


MAGNUS  FORCE  PRODUCTION  ON  SLENDER  BODIES  OF  REVOLUTION 

AT  LOW  ANGLES  OF  ATTACK 


To  study  three-dimensional  flow  about  a  spinning  body 
of  revolution  the  full  viscous  equations  of  motion  must  be 
solved.   The  general  equations  of  motion  predict  that  the 
side  force  is  a  function  of  angle  of  attack,  non-dimensional 
spin  rate  (the  ratio  of  the  surface  speed  CO  r  and  the  free- 
stream  velocity  U) ,  body  fineness  ratio,  Reynolds  number, 
and  Mach  number.   The  complete  equations  of  motion  are  much 
too  difficult  to  solve,  hence  simplifications  must  be  made. 
The  simplest  problem  involves  incompressible  flow  with  a 
laminar  boundary  layer,  a  small  angle  of  attack,  and  a  small 
non-dimensional  spin  rate. 

A  useful  method  for  studying  such  a  flow  is  by  using  a 
cross-flow  analogy.   The  basic  assumption  made  to  simplify 
the  problem  is  that  the  flow  along  the  body  is  independent 
of  the  flow  across  the  body.   This  assumption  is  most 
likely  inaccurate  for  high  angles  of  attack  and  low  body 
fineness  ratios.   Experimental  test  data,  however,  should 
define  the  region  of  validity  for  the  application  of  this 
theory. 

Using  the  cross-flow  analogy,  a  spinning  missile  placed 
in  a  flow  at  an  angle  of  attack  is  assumed  to  have  the  flow 
velocity  about  it  resolved  into  two  components  as  shown  in 
Figure  1. 
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Figure  1.   Velocity  Components  over  a  Spinning 
Missile 


The  drift  experienced  by  the  missile  during  flight  can 
be  explained  by  investigating  the  effect  of  the  cross-flow 
component  (the  velocity  component  perpendicular  to  the 
missile) . 

From  inviscid  flow  theory  it  is  known  that  lift  per 
unit  span  on  a  spinning  cylinder  in  a  uniform  stream  is 

L   =  pUT   . 
where  T   is  the  circulation  of  the  flow.   The  flow  pattern 
for  the  spinning  cylinder  in  inviscid  flow  is  shown  in 
Figure  2. 


Figure  2.   Streamline  Pattern  of  a  Spinning 
Cylinder  in  Inviscid  Flow 
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Inviscid  flow  theory  is  helpful  for  understanding  how 
lift  is  generated  on  a  body  in  a  flow,  but  it  ignores  how 
the  circulation  responsible  for  the  lift  is  generated. 

Real  flow  is  viscous  and  does  not  behave  as  described 
in  inviscid  flow  theory  near  the  body  surface.   From  bound- 
ary layer  theory  it  is  known  that  the  flow  around  a  spinning 
cylinder  forms  a  boundary  layer.   This  boundary  layer  is 
ultimately  shed  somewhere  on  the  rear  side  of  the  cylinder 
in  the  flow.   Because  of  the  cylinder's  rotation,  the 
boundary  layer  formed  is  not  symmetric  about  the  cylinder 
nor  is  it  shed  symmetrically  in  the  flow.   The  boundary  layer 
is  shed  in  the  form  of  vortices  and  is  shown  schematically 
in  Figure  3. 


"r  =  Ua 


Figure  3.   Boundary  Layer  and  Shed  Vorticity  for 
a  Spinning  Cylinder  in  Inviscid  Flow 


The  net  vorticity  is  in  the  direction  opposite  to  the 
rotation  of  the  cylinder,  in  this  case  counter-clockwise. 
Because  the  total  vorticity  must  be  zero  for  the  flow  as 
a  whole,  a  circulation  in  the  clockwise  direction  is  gen- 
erated and  it  is  this  circulation  which  provides  a 
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major  contribution  to  the  Magnus  force  acting  on  the 
cylinder. 

In  addition  to  the  circulation,  boundary  layer  dis- 
placement thickness,  radial  pressure  gradient  in  the  bound- 
ary layer,  and  the  presence  of  skin  friction  shearing 
stresses  contribute  to  the  production  of  a  Magnus  force. 
The  displacement  thickness  is  the  distance  by  which  the 
boundary  layer  would  have  to  be  displaced  if  the  entire 
flow  were  imagined  to  be  f rictionless .   A  representation 
of  the  displacement  thickness  is  shown  in  Figure  3.   It 
is  not  symmetric  about  the  cylinder  because  of  the  effect 
the  cylinder's  spin  has  on  the  boundary  layer  growth. 
Knowledge  of  the  displacement  thickness  permits  use  of  a 
"displaced"  body  in  place  of  the  actual  body  such  that  the 
flow  around  the  body  can  be  considered  f rictionless .   As 
can  be  seen  in  Figure  3,  the  magnitudes  of  the  velocity 
components  at  the  boundary  layer  edge  of  the  advancing  sid 
of  the  cylinder  are  less  than  those  of  the  retreating  side 
and  hence  the  outer  inviscid  flow  velocity  is  less  below 
the  cylinder  than  above.   Because  this  outer  flow  is 
frictionless ,  Bernoulli's  equation  predicts  that  the 
pressure  below  the  cylinder  is  less  than  that  above,  and 
a  lift  is  generated. 

Although  small,  it  can  be  shown  that  there  is  a  radial 
pressure  gradient  in  the  boundary  layer  around  the  cylinder 
due  to  the  surface  curvature.   When  integrated,  the  change 
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in  surface  pressure  due  to  this  pressure  gradient  contri-* 
butes  to  the  cylinder's  lift. 

The  presence  of  friction  on  the  surface  of  the  cylinder 
is  also  a  cause  of  lift.   When  the  frictional  shearing 
stress  is  integrated  over  the  cylinder's  surface  there  is 
a  component  of  the  resultant  force  in  the  direction  of 
lift. 

Although  circulation  is  a  major  factor  in  creating 
lift  on  the  spinning  cylinder,  the  effects  of  the  radial 
pressure  gradient  and  surface  friction,  however  small, 
should  be  included  in  the  analysis  of  the  total  Magnus 
force  developed  by  the  cylinder. 
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DESCRIPTION  OF  APPARATUS 

Testing  to  obtain  experimental  force  data  was  done  in 
the  Naval  Postgraduate  School  low  speed  wind  tunnel.   The 
tunnel  test  section  was  hexagonal  in  shape  with  maximum 
dimensions  of  3.5  x  5.0  feet. 

The  model  used  for  testing  was  a  five-sectioned 
aluminum  cylinder  2.950  inches  in  diameter  with  an  ogive 
nose  shape  as  shown  in  Figure  4.   The  missile  was  con- 
structed in  sections  so  that  length-to-diameter  ratios  of 
9.9,  7.8,  and  6.0  could  be  tested. 

Figure  5  is  a  line  drawing  of  the  model  and  sting 
balance  system.   The  missile  was  supported  by  two  Fafnir 
KD66  precision  ball  bearings  which  were  press-fitted  onto 
the  sting  balance.   Each  bearing  was  enclosed  in  a  stainless 
steel  collar  which  allowed  the  missile  to  be  fastened  to 
it  with  screws.   Mounted  at  the  end  of  the  balance  was  a 
1/8  horsepower  Standard  Pneumatic  rotary-vane,  bi-directional 
compressed  air  motor.   Figure  6  shows  the  balance,  bearing, 
and  motor  configuration.   The  compressed  air  was  supplied 
to  the  motor  through  1/4-inch  copper  tubing  routed  along 
the  missile  support  apparatus. 

Mounted  on  the  drive  shaft  of  the  air  motor  was  an 
aluminum  drive  disk.   Located  near  the  disk  was  an  elec- 
tromagnetic tachometer  pickup  which  sensed  passage  of 
ferrous  material  (steel  screws)  in  the  aluminum  disk. 


15 


u 


o 

<D 

00 

< 
> 

•H 
•M 
•H 

o 
+-» 

03 


CD 
O 

S 

<D 


DO 

•H 


16 


U 

e 

PQ 


P. 
co 

•  H 
< 


U 

c 

as 
PQ 

C 


o 

O 

s 

•H 


o 

m 
o 

•H 

<D 
Q 


<D 

•H 


17 


O 
+-> 

o 

s 

u 

•H 

< 

C 

cd 

•H 
*H 

Ctf 


u 

c 

i-H 


bO 
•H 

PL, 


18 


Also  mounted  on  the  balance  were  two  four-element 
strain  gage  bridge  circuits  which  were  used  to  measure 
strain  in  the  sting  balance  and  thus  provide  the  means  of 
determining  the  Magnus  force  and  its  center  of  pressure. 
Each  gage  was  a  Baldwin-Lima-Hamilton  strain  gage,  model 
number  FAE  25  12S13,  with  a  gage  factor  of  2.08.   The 
arrangement  of  the  gages  for  one  bridge  is  shown  in 
Figure  7.   The  gages  were  wired  in  such  a  manner  as  to 
sense  only  side  force  and  the  Magnus  moment  applied  to  the 
balance.   Because  of  the  beam  geometry  and  the  longitudinal 
arrangement  of  the  gages  there  was  no  response  to  any 
applied  torque.   Lift  forces  on  the  beam  also  had  no  effect 
on  the  bridge  output  because  of  gage  arrangement.   The 
bridge  remained  balanced  when  lift  forces  were  applied 
because  one  pair  of  gages  (#2  and  #4  as  shown)  was  in  com- 
pression while  the  other  pair  (#1  and  #3)  was  in  tension. 
The  currents  in  each  pair  were  different  because  resist- 
ances were  different,  but  because  the  voltage  drop  was  the 
same  across  each  pair  the  bridge  remained  balanced. 
Because  each  bridge  circuit  was  made  up  of  four  active 
elements  the  strain  indication  was  magnified  by  a  factor 
of  four.   The  bridge  was  also  temperature  compensating 
because  all  elements  were  affected  the  same  by  changes  in 
temperature.   Figure  8  is  a  picture  of  one  side  of  the 
balance  with  mounted  strain  gages. 

The  balance  was  screwed  at  its  base  into  a  tapered, 
circular  cross-section,  stainless  steel  sting.   The  sting 
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was  held  in  place  at  its  base  by  a  support  and  could  be 
rotated  90°  for  measurements  in  either  the  horizontal  or 
vertical  plane.   A  90  x  6  x  1-inch  aluminum  sting  support 
was  bolted  top  and  bottom  to  the  wind  tunnel  to  support 
the  missile  and  sting  arrangement.   The  beam  Avas  machined 
so  that  the  entire  sting  could  be  rotated  in  the  vertical 
plane  to  provide  angle  of  attack  as  shown  in  Figure  4. 
The  sting  support  was  raised  or  lowered  to  keep  the  model 
centered  in  the  test  section.   The  angle  of  attack  was 
measured  with  a  clinometer. 

Figure  9  is  a  wiring  diagram  of  the  electronics  used 
to  obtain  data.   An  input  voltage  of  6  volts  was  supplied 
to  each  bridge  circuit  by  a  Moxon  Electronics  SRC  Division 
power  supply,  model  3564.   Each  bridge  output  was  amplified 
by  a  factor  of  100  by  an  Aztec  model  886  DC  amplifier  and 
recorded  on  a  Hewlett-Packard  two  channel  chart  recorder, 
model  7100B.   Two  potentiometers  per  bridge  circuit  were 
used  to  initially  balance  each  bridge  before  testing  was 
begun.   Figure  10  is  a  picture  of  the  apparatus  used  in 
the  testing.   In  addition  to  the  digital  counter,  a 
General  Radio  strobe,  model  1540-P1,  was  also  used  to 
determine  missile  spin  rate  to  the  nearest  10  RPM. 
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CALIBRATION  AND  DATA  REDUCTION  METHODS 

Airspeed  in  the  tunnel  test  section  was  measured  using 
water  manometry.   A  pitot  tube  which  measured  total  and 
static  pressures  in  the  test  section  provided  the  inputs 
to  the  water  manometry  set  up.   Because  the  pitot  tube  was 
located  in  the  most  forward  part  of  the  test  section  and 
only  six  inches  off  the  test  section's  floor,  calibration 
of  the  system  was  necessary.   To  calibrate  the  tunnel's 
system  another  pitot  tube  was  mounted  in  the  center  of  the 
test  section  and  readings  were  compared  throughout  the 
range  of  airspeeds  possible  in  the  tunnel.   These  corrections 
to  tunnel  readings  were  made  in  the  computer  data  reduction 
program. 

In  addition  to  pitot  tube  calibration,  corrections  to 
the  test  section  velocities  were  necessary  due  to  test  sec- 
tion blockage  by  the  model.   The  correction  applied  is 
developed  in  Pankhurst  and  Holder   [9] : 

V+     =  V  ,  (1.0  +  c  ) 
true    obs v       s' 

where  V  »   is  the  velocity  directly  measured  in  the  tunnel 
test  section,  €   is  the  tunnel  blockage  factor  and  is  given 


by: 


IT   h  Vol 

€c  =  C-4")2  T  (1.0  +  .4D/L) 


c3/2 

IT  U 

where  (-j-)  2  is  a  factor  appropriate  to  a  body  of  revolu- 
tion, T   is  a  factor  dependent  on  the  shape  of  the  tunnel 
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cross  section,  Vol  is  the  volume  of  the  model,  D/L  is  the 
reciprocal  of  the  model  fineness  ratio,  and  C  is  the  tunnel 
test  section  cross  sectional  area. 

Corrections  to  test  data  due  to  the  effect  of  air  hose 
pressurization  on  bridge  circuit  output  were  necessary. 
Pressure  applied  to  the  air  motor  stiffened  the  rubber  hoses 
leading  to  the  motor  and  caused  a  deflection  in  the  sting 
balance  which  in  turn  affected  strain  gage  bridge  output. 
Corrections  for  pressurization  were  made  in  the  computer 
data  reduction  program  to  correct  for  these  effects.   Input 
data  for  the  pressurization  corrections  were  obtained  by 
prior  calibration  of  the  sting  balance  which  gave  chartline 
deflection  for  any  given  air  pressure. 

The  sting  balance  system  described  previously  was  cali- 
brated by  hanging  a  known  weight  a  known  distance  from  each 
bridge  circuit  location  and  recording  each  bridge's  output 
in  units  of  chart  lines  on  the  Hewlett-Packard  chart 
recorder.   Calibration  apparatus  is  shown  in  Figure  11. 
By  plotting  the  moment  versus  chartline  deflection  for 
each  bridge  circuit  linear  graphs  were  obtained  whose  slopes 
were  calculated.   Test  results  (chartline  deflections  for 
unknown  forces)  were  multiplied  by  the  slopes  for  each 
circuit  to  determine  moments.   By  knowing   the  moment  at 
each  bridge  circuit  the  applied  force  and  its  location  as 
shown  in  Figure  12  were  calculated  using  the  following 
simultaneous  algebraic  equations: 
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Figure  12.   Strain  Gage  Bridge/Beam  Geometry 


F  =(M   -  M2)/d 


X   =  M0/F 
cp    2 
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TESTING  PHILOSOPHY  AND  PROCEDURE 

Testing  to  determine  the  Magnus  moments  was  done  with 
the  objective  of  obtaining  data  to  verify  the  theoretical 
equation  developed  by  Power  in  1971   [8] ,  which  predicts 
the  Magnus  force  coefficient  developed  by  a  spinning  body 
of  revolution.   Power's  equation  showed  that  the  Magnus 
force  coefficient  was  dependent  on  angle  of  attack,  non- 
dimensional  spin  rate  (ratio  of  Cylinder  surface  speed  to 
flow  velocity),  model  length-to-diameter  ratio,  the 
Reynolds  number  based  on  missile  length,  and  the  cross  flow 
Reynolds  number. 

1   "   v  *L 


Cy  =<|ks  C^)2  ^60.9757  +  ^h    (jp  [256  C±   R 


357.37]  -  227-0112  (L)2\  (l) 

R       D   J  K    J 

Lt 


.6q  f 
where   C,  =  . rr   _  and  f,  given  in  Figure  12,  is  a  functio 

of  cylinder  spin  rate  and  angle  of  attack.   By  varying 

angle  of  attack,  non-dimensional  spin  rate,  and  the  model 

length-to-diameter  ratio  while  maintaining  constant  Reynolds 

numbers  for  each  angle  of  attack,  data  could  be  obtained 

for  three  length-to-diameter  ratios  and  various  spin  rates 

which  would  enable  the  theoretical  equation  to  be  compared 

with  experimental  data.   Because  the  non-dimensional  spin 

rate  was  in  part  dependent  upon  the  Reynolds  number  based 
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on  model  length,  two  Reynolds  numbers  were  chosen  for  testing 
which  were  compatible  with  all  three  model  lengths  and  wind 
tunnel  capabilities.   The  two  Reynolds  numbers  selected  were 
700,000  and  1,600,000. 

It  was  found  that  pre-determining  non-dimensional  spin 
rate  (p)  for  the  given  Reynolds  number  was  impractical 
because  the  spin  rate  required  for  each  p  was  difficult  and 
time-consuming  to  set  during  testing.   Therefore,  the  pro- 
cedure decided  upon  was  to  pre-select  various  air  pressures 
to  the  air  motor  and  accept  whatever  spin  rates  and  hence 
p's  that  resulted. 

Procedure  played  a  major  part  in  the  testing  because 
of  the  time  that  could  be  saved  through  proper  planning. 
Because  two  tunnel  velocities  and  a  wide  range  of  non- 
dimensional  spin  rates  (plus  and  minus)  were  tested  for 
each  angle  of  attack,  testing  was  extremely  time  consuming. 
By  establishing  a  procedure  which  involved  a  minimum  number 
of  configuration,  tunnel  velocity,  and  model  spin  rate 
changes,  the  time  consumed  for  any  part  of  the  testing  was 
reduced  to  a  minimum. 

Upon  completion  of  testing  all  raw  data  obtained  were 
typed  on  computer  cards  and  fed  into  an  IBM  360  digital 
computer  and  reduced  by  the  data  reduction  program. 
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DATA  ANALYSIS 

Due  to  various  physical  phenomena  encountered  during 
the  testing  the  raw  data  obtained  were  very  difficult  to 
analyze  to  obtain  meaningful  results.   A  sample  of  the 
type  of  data  obtained  during  the  testing  from  the  Hewlett- 
Packard  chart  recorder  is  shown  in  Figure  14. 

The  raw  data  indicated  that  the  missile  vibrated  in 
two  modes  during  the  testing.   The  amplitude  of  the  short 
period  mode  varied  between  three  and  five  chartlines 
(three  and  five  millivolts)  and  appeared  to  be  at  a  con- 
stant frequency  throughout  the  testing.   The  amplitude  of 
the  long  period  mode  varied  between  eight  and  eleven  chart- 
lines  in  amplitude  (measured  from  the  mean  of  the  short 
period  oscillation)  but  its  frequency  was  random. 

Data  were  obtained  by  considering  the  long  period  mode 
first.   The  best  estimate  of  its  average  amplitude  was 
made  and  this  information  was  modified,  if  necessary,  by 
considering  the  superimposed  average  short  period  amplitude 
In  this  way  a  constant  voltage  value  was  determined  for  use 
in  the  data  reduction  program.   Because  of  the  constant 
frequency  of  oscillation  of  the  short  period  mode,  averag- 
ing was  easy  compared  to  averaging  the  random  oscillations 
of  the  long  period  mode.   These  vibrational  modes  were 
caused  by  at  least  four  major  effects.   First,  tunnel  wall 
vibration  was  transmitted  to  the  model  through  the  model 
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Figure  14.   Sample  of  Raw  Data  Output 
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support  system.   The  sting  and  support  were  bolted  directly 
to  the  tunnel  floor  and  ceiling  in  the  test  section.   Lack 
of  time  prevented  the  construction  of  the  originally 
designed  isolated  support  system.   Second,  mass  unbalance 
of  the  missile  itself  caused  a  vibration  when  the  missile 
was  rotating.   Even  though  the  missile  and  bearing  collars 
were  finely  machined  there  was  a  mass  unbalance  due  to 
varying  missile  wall  thickness  and  bearing  slop.   Third, 
asymmetric  vortex  shedding  and  leeside  separation  of  the 
flow  about  the  missile  added  to  the  vibration.   This  was  a 
function  of  Reynolds  number  and  rotation  rate.   Fourth,  the 
tunnel  test  section  velocity  was  not  constant  and  large  air 
flow  fluctuations  occurred.   The  solution  to  this  problem 
would  be  to  install  flow-straightening  screens  to  the  tunnel 
to  ensure  constant  low  turbulence  flow  in  the  test  section. 

During  the  testing  a  breakdown  forced  the  removal  of 
one  of  the  two  counter-rotating  Bendix  constant  speed  pro- 
pellers used  to  generate  the  tunnel's  velocity.   This 
severely  reduced  the  maximum  velocity  that  could  be 
attained  and  subsequently  delayed  the  testing  programs  until 
maintenance  was  completed  and  a  new  set  of  parameters  for 
testing,  reflecting  the  change  in  tunnel  speed  capability, 
could  be  computed.   The  removal  of  the  propeller  was  nec- 
essary because  of  lack  of  parts.   The  long  period  of  vibra- 
tion was  attributed  primarily  to  the  tunnel  test  section 
velocity  fluctuations  while  the  short  period  oscillations 
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were  attributed  to  tunnel  wall  vibration,  mass  unbalance 
in  the  missile,  asymmetric  vortex  shedding,  and  flow 
separation. 

Because  of  the  magnitude  of  the  oscillations  and 
tunnel  flow  variations,  error  bounds  were  computed.   Due  to 
the  type  of  strip  chart  paper  used  (ten  chart  lines  per 
inch)  it  was  assumed  that  the  raw  data  were  accurate  to 
within  ±1  chartline.   It  was  also  estimated  that  tunnel 
flow  fluctuation  was  ±1  centimeter  of  water  at  dynamic 
pressures  above  10  centimeters  of  water  and  ±.5  centimeter 
below  dynamic  pressures  of  10  centimeters.   Evaluation  of 
the  error  bound  for  each  data  point  was  done  by  computer 
analysis  and  representative  values  are  shown  in  Figure  15. 
The  error  in  C   due  to  tunnel  velocity  fluctuations  was 
about  11%  at  each  Reynolds  number  tested.   The  major  error 
contribution  was  due  to  the  vibration  and  the  resulting 
inability  to  interpret  the  strain  gage  readings  accurately. 
The  uncertainty  bound  for  each  data  point  was  computed  by 
adding  one  chartline  to  the  aft  bridge  circuit  reading  and 
subtracting  one  chartline  from  the  forward  bridge  circuit 
reading.   In  this  way  the  maximum  testing  error  could  be 
calculated.   It  can  be  seen  in  Figure  15  that  the  larger 
of  these  uncertainties  occurred  for  the  lowest  Reynolds 
number  and  highest  angles  of  attack  (C   is  directly  pro- 
portional to  angle  of  attack  and  inversely  proportional 
to  tunnel  velocity) .   Because  of  the  magnitude  of  the  error 
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Figure  15.   Maximum  Error  Bound  Samples 
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bounds,  interpretation  of  the  data  was  very  difficult. 
Similar  problems  have  plagued  all  experimenters  investigat 
ing  Magnus  forces  to  date.   The  forces  measured  are  very 
small  and  the  data  accuracy  required  is  better  than  that 
usually  required  for  measuring  drag  forces. 
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RESULTS 

A  representative  summary  of  the  data  obtained  during 
experimentation  is  presented  in  tabular  form  in  Appendix  B 
and  graphically  in  Figures  16  through  27.   Figures  16 
through  20  are  plots  of  corrected  data.   At  each  angle  of 
attack  the  data  were  reduced  by  subtracting  the  Magnus 
force  coefficient  obtained  for  the  no-spin  condition  from 
each  data  point.   This  ensured  that  the  curve  of  force 
coefficient  versus  spin  rate  would  pass  through  the  origin 
as  predicted  by  theory.   This  procedure  corrected  the  data 
for  the  experimental  side  force  coefficient  resulting  from 
any  missile  yawing  mis-alignment  in  the  tunnel.   Represen- 
tative error  bounds  for  the  data  points  are  shown  for  each 
of  the  first  five  plots  (16  through  20) .   As  discussed 
earlier,  the  bound  is  larger  for  the  smaller  Reynolds 
number  and  higher  angles  of  attack.   Small  Q  and  p  theory 
predicts  that  a  plot  of  Cy   versus  p  is  linear.   Hence  the 
corrected  data  points  in  Figures  16  through  20  were  fitted 
with  a  straight  line  whose  slope  d^y/dp   was  cross-plotted 
versus  angle  of  attack  for  each  length-to-diameter  ratio 
and  Reynolds  number.   This  slope  was  then  adjusted  to 
yield  a  smooth  variation  for  dCy/dp  versus  Ct  .  Justification 
for  this  method  was  provided  by  the  size  of  the  error  bound 
in  each  case. 
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Figures  21  through  25  are  the  resulting  plots  of  the 
Magnus  force  coefficient  Cy  versus  angle  of  attack  Q  .   Also 
shown  are  results  developed  by  Power  and  Kelly  for  calcu- 
lating Magnus  force  coefficients.   Figures  26  and  27  are 
plots  of  Magnus  force  coefficient  versus  length-to-diameter 
ratios  for  constant   Reynolds  number  and  angle  of  attack 
and  for  various  spin  rates. 
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Figure  21. 


Magnus  Force  Coefficient 
vs.  Angle  of  Attack 
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Figure  22. 


Magnus  Force  Coefficient 
vs.  Angle  of  Attack 


46 


.7   r  C 


p  =  .4 


P  =  -2 


=  .2 


Figure  23. 


Magnus  Force  Coefficient 
vs.  Angle  of  Attack 
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Figure  24. 


Magnus  Force  Coefficient 
vs.  Angle  of  Attack 
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Magnus  Force  Coefficient 
vs.  Angle  of  Attack 
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Figure  26.   Magnus  Force  Coefficient  vs.  L/D 
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CONCLUSIONS 

The  graphical  data  suggest,  but  with  the  error  bounds 
of  the  experimentation  cannot  verify,  that  the  Magnus 
force  might  be  linear  in  both  Ct    and  p  for  small  angles  of 
attack  and  small  spin  rates  at  constant  Reynolds  numbers. 
Other  experimenters  have  found  non-linear  variations  of  Cy 
at  angles  of  attack  greater  than  10°.   Power's  theory  pre- 
dicts a  linear  variation  with  p  and  non-linear  variation 
with  CL     while  Kelly's   [7]  predicts  a  linear  variation  with 
both  p  and  a  .   However,  due  to  the  magnitude  of  the  error 
associated  with  each  data  point  obtained,  a  great  deal  of 
latitude  for  data  interpretation  was  available  once  the 
data  were  reduced.   The  data  scatter  indicates  that  if  the 
sources  of  error  discussed  in  the  data  analysis  section 
could  be  reduced,  the  results  would  be  more  valid  for 
comparison. 

In  general  the  forces  obtained  were  much  lower  than 
those  predicted  by  Kelly's  or  Power's  theory  as  shown  by 
Figures  21  through  25.   These  theories  were  developed 
assuming  an  open-ended  cylinder  while  the  model  had  an 
ogive  nose  shape.   This  difference  alone,  however,  cannot 
account  for  the  discrepancies  between  theory  and  experi- 
mental results.   Figures  26  and  27  indicate  that  as  the 
length-to-diameter  ratio  of  the  missile  increases,  so  does 
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the  side  force.   However,  comparison  of  the  experimental 
data  with  theory  shows  that  the  variation  of  Cy  with  L/D 
in  the  region  of  these  experiments  is  smaller  than  predicted 
by  theory. 

The  experimental  data  presented  in  Figures  21  through 
25  were  used  to  calculate  C,  in  equation  (1).   The  results 
showed  C,  to  be  a  small  negative  number.   This  implies 
that  flow  circulation  contributed  negatively  to  the  Magnus 
force,  which  contradicts  physical  reasoning.   If  C,=0,  the 
resulting  equation  is  essentially  what  Kelly  predicted. 
However,  the  theory  must  not  be  discounted  simply  because 
the  experimental  data  do  not  support  it.   When  the  error 
associated  with  the  data  is  considered  it  is  readily 
apparent  that  the  data  may  not  correlate  with  theory. 

Because  the  measurement  of  these  forces  is  so  critical, 
a  better  method  for  measuring  them  must  be  devised. 
Although  the  method  used  did  provide  some  meaningful  data 
which  supported  theory  to  some  extent,  the  necessity  for 
finer  instrumentation  and  refinement  of  equipment  is 
obvious.   Precise  machining  was  necessary  for  all  of  the 
apparatus,  but  no  amount  of  machining  or  expensive,  elab- 
orate electronic  equipment  will  enable  accurate  data  to 
be  obtained  if  other  parameters  such  as  poor  tunnel  flow 
characteristics  are  a  part  of  the  system.   The  model,  as 
constructed,  can  be  tested  with  a  variety  of  nose  shapes 
so  that  comparisons  can  be  made  of  the  effects  different 
forebody  shapes  have  on  the  production  of  Magnus  forces. 
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However,  until  more  accurate  results  are  obtainable  it 
would  be  useless  to  test  different  forebody  shapes.   Be- 
cause the  data  obtained  do  indicate  that  some  side  force 
was  measured  there  is  sufficient  reason  for  continuing  the 
experimentation  with  re-designed  equipment  to  obtain  more 
meaningful  results. 
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RECOMMENDATIONS  FOR  FURTHER  STUDY 

i 

1.  Design  a  more  sensitive  balance  that  is  less 

affected  by  air  motor  hose  pressurization. 

2.  Redesign  the  sting  mounting  so  that  it  is  iso- 
lated from  test  section  vibrations. 

3.  Install  flow  straightening  screens  to  reduce 
tunnel  test  section  flow  fluctuations. 

4.  Re-examine  the  circulation  coefficient  C   to 
correlate  theory  by  taking  new  data  with  improved  exper- 
imental system. 

5.  Perform  nose  shape  studies  to  assess  the  affects 
of  forebody  shape  on  Magnus  force. 
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COMPUTFR  PROGRAM  TO  REDUCE  RAW  DATA  TO  DETERMINE  MAGNUS 
FORCE  AND  COEFFICIENT  WITH  ERROR  BOUNDS 

DIMENSION  ALISTt 100),BLIST( 100) .FLIST(IOO) tGLIST( 100) , 
PLISTUOO) ,CLIST(100) 


EXPLANATION  OF 
XL/XLF   - 
DF 
XLD 
T/TR 
P 
RHO 
XMU 
R 
G 
XMF/XMA   - 

DF 
AOA/ARAD  - 
0 
V/VTRU  - 

EPS 

TSEC 

XCO 

CONE 

VOL 

REV/XO 

REL/REC   - 

XCO 
S 

FCOR/ACOR  - 


SYMBO 
MISS 
MI  SS 
MISS 
TEMP 
AT  MO 
ATMO 
COEF 
GAS 
GRAV 
SLOP 
FROM 
DI  AM 
ANGL 
DYNA 
TUiMN 
VELO 
TEST 
LENG 
VOLU 
VOLU 
REVO 
SECO 
REYN 
FLOW 
LENG 
MISS 
CORR 
HOSE 


LS  U 
ILE 

I  LE 

ILE 

ERAT 

SPHE 

SPHE 

FICI 

CONS 

ITY 

ES  C 

STI 
ETER 
E  OF 
MIC 
EL  V 
CITY 

SEC 
TH  0 
ME  0 
ME  0 
LUTI 
NO 
OLDS 


SED  I 
LENGT 
DI  AME 
LENGT 
URE  ( 
RIG  P 
RIC  D 
ENT  0 
TANT 
CONST 
F  MOM 
NG  CA 
OF  M 
ATTA 
PRESS 
ELOCI 
CORR 
TION 
F  MIS 
F  MIS 
F  MIS 
ONS  P 


N  THE  PROGRAM 
H  ( IN./ET. ) 

TER  IN  FEET 

H-TO- DIAMETER  RATIO 

FAHF  ENHEIT/RANKINE) 

RESSURE  (IN.  HG.) 

ENS  ITY 

F  VISCOSITY  (AIR) 


ANT 

ENT 

LIB 

ISS 

CK 

URE 

TY 

ECT 

CRC 

SIL 

SI  L 

SIL 

ER 


VS  CHART  DEFLECTION 
RATION 
ILE  (FT.) 
(DEGREES/RADIANS) 

IN  CM  H20 
( MEASURED/CORRECTED) 
I  ON  FACTOR 
SS  SECTIONAL  AREA 
E  NOSE  SHAPE 
E  NOSE  PIECE 

c 

MINUTE/RADIANS  PER 


NUMBERS  BASED  ON  LENGTH/CROSS 


TH  OF  MISSILE  NOSE  SHAPE 
ILE  CROSS  SECTIONAL  AREA 
EC  71 ONS  TO  DATA  DUE  TO  EFFECTS  OF 
STIFFENING  ON  STING 


WRITE (6,  10) 


READ    IN    CIRCULATION    FUNCTION    DATA    FOR 
THEORETICAL    EQUATION     FOR    MAGNUS     FORCE 


USE     IN    POWER'S 
COEFFICIENT 


RE  AD (5,  14) (PL  I  ST (J )  ,J  =  1,17) 
READ(5,14)(CLIST(J),J=1,17) 


DO    100    1=1, 10 


READ(5.1 

IF(E.EO. 

IF(E.GT. 

WRITE(6, 

GO    TO    2 

WRITF(6, 


1 )L,NDATE,M,E 
1.0)G3    TO    9 
5.0)G3     TO     1 
12JNDATE 

13JNDATE 


RFAD  IN  DATA  FOR 
AIR  PRESSURES 


FORCE  READING  CORRECTIONS  DUE  TO  HOSE 


READ(5,14)(ALIST(J)  ,J=1,L) 
RFAD(5,14)(BLIST(J),J=1,L) 
R^ AD (5. 14)  (FLIST (J  )  ,J=1,L> 
READ( 5, 14) (GLI  ST( J)  ,J=1, L) 


DO  101  N=1.M 

READ  IN  MISSILE  LENGTH, 
FAHRFNHEIT),  ATMOSPHERIC 
PARAMETER  USED  TO  DET~ 
PRESSURE 


TEMPERATURE  (DEGREES 
PRESSURE,  AND  A  MISSILE 
ERMINE  LOCATION  OF  CENTER  OF 


READ(5,  15) XL 
XLF  =  XL/12 


T,P,D,K 
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DF    =    2.950/12. 
XLD    =     XLF/DF 
WRITE(6. 16 )XL, XLD 

TR    =    T+459.67 
R    =    53.34 
c    "~     -^  'y    t  ~7  / 

RHG    =    P*2116.22/(29.92*G*R*TR) 

X    =    2.26967    E-08 
XMU    =    X=:  (TR-: -M..5)  /(TR+198.72) 
XMF    =    63.7/274.8 
XMA    =    86.4/323.7 

DO    102     J=1,K 

RFAD    IN    MISSILE     ANGLE    OF    ATTACK,     DYNAMIC     PRESSURE    IN 
CM    H?0    (CORRECTED    TO    TRUE    0    IN    PROGRAM),     THE 
REVOLUTIONS    PER    MINUTE,     AND    THE    TWO    CHARTLINE 
READINGS     (1    CHARTLINE    PER    MILLIVOLT) 

READ(5,17) AOA, 0, REV , VF , VA, PHOS E 
SPIN    =     REV 
ANGL    =    ABS(AOA) 
ARAD    =    ANGL'  3. 14159/180. 
CA    =    COS(ARAD) 
SA    =    SIN(ARAD) 
Z    =     0*62. 4-'- 1  .133/  (RHO- 15.24) 
V    =     SORT(Z) 

INSFRT    DATA    TO    DETERMINE    VELOCITY    ERROR    BOUNC 

IF(O.GT.IO.OO) GO    TO    3 

OX    =    0+0.5 
GO    TO    4 

3  OX    =    0+1.0 

4  ZP    =    0X-:<62.4M.133/(RH0<  15.24) 
VP    =    SORT(ZP) 

14.41 

.01973 
7.244/12. 

CONF+3. 14159*(DF**2)* (XLF-XCO) /4. 
.75*>(  (3. 14 159/4.)**. 5)1  (  1  .+.4*DF/XLF  J*  VOL/ 
**1.5) 
V-/(  l.+EPS) 
VP>:  {  l.+  EPS) 
REVr-2.  £3. 14159/60. 
DF*XO/VTRU 
DF*  XC/VTRP 
RH03  XLF*  VTRU-'-CA/XMU 
R  HO -••  D  F>  V TR  U>  SA  /  X M U 
(  (DF/2.  )•  -2):-3.14159 

PIF2(PH0SE, ALIST,L,BLIST) 

VF-FCOR 

PIF2(PH0SE,FLI ST,L,GLIST) 

VA-ACCR 

XMF*  VFC 

XMA-rVAC 

DATA    INPUT    FOR    DETERMINING    ERROR    IN    INITIAL    RAW    DATA 
READINGS 


TSEC  = 

CONE  = 

XCO  = 

VOL  = 

EPS  = 

(TSEC 

VTRU  = 

VTRP  = 

5 

XO  = 

6 

PB  = 

PBP  = 

REL  = 

REC  = 

S  = 

FCOR  = 

VFC  = 

ACOR  = 

VAC  = 

XMOF  = 

XMCA  = 

VFP 

= 

VFC+1.0 

VFM 

= 

VFC- 1.0 

VAP 

= 

VAC  +  1  .0 

VAM 

= 

VAC-1.0 

XMFP 

= 

XMF" -VFP 

XMFM 

= 

XMF=:  VFM 
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XMAP    = 
XMAM    = 


XMA^VAP 
XMA*VAM 


CALCULATION    OF    MAGNUS 
ERROR    BOUND 


FORCE,     FORCE    COEFFICIENT,     AND 


Y    =  (XMOA-XMCF) /5 . 

XCP    =  XMOF/Y 

CY    =  2.*Y/(RH0--VTRU?.*VTRU*'S) 

YP    =  ( XMAM-XMFP) /5. 

CYP    =  2.*YP/{  RHC'VTRU^VTRU^-S) 

YM    =  (XMAP-XMFMJ/5. 

CYM    =  2  .•'•  YM/  (  RH  0*  V  TRU-  VTRU*  S ) 

CYVP    =  2.-':Y/(RHO'-VTRP*  VTRP'.:  S) 


IF(A0A.LT.1.0)G0    TO    7 
IF ( PR.  FO. 0.0) GO    TO    7 
PSIN    =    ABS(PB/SA) 

PIF2(PSIN,PLIST,17,CLIST) 

.6-"  APAD-'-F/(  128. *3. 14159*  PB) 

(  ARADO  PB'XLD"  XL  D/  (R  EL**  .5  )  )$(  60  .97  +  (  256.1-C  1# 

.  7  5 )  -  3  5  7  .  3  7 )  ■-> -XL  D  /  (  RE  L*  * .  5 )  -2  2  7 .  0  lr  X  LD*  XLD/  RE  L  ) 


F    = 
CI    = 
CYPW    = 
(REC** 

GO    TO    8 
CYPW    =    0 


0 


PRINTOUT    DATA    DESIRED 

8    WRITE (6. 18)A0A,REL, VTRU, SPIN,PB,Y,CY,CYP,CYM,CYVP,CYPW 

10?  CONTINUF 
101  CO NT  IN UE 
100    CONTINUF 

FORMAT    CONTROL    STATEMENTS    ARE: 


10 

11 
12 
13 
14 
15 
16 


17 
18 

19 


Bf ///T31,' MAGNUS     FORCE    AND1 


1  ) 

•TESTING 
•TESTING 


DATE 
DATE 


I  3, IX, ' JANUARY  1973') 
13, IX, 'FEBRUARY  1973' 


FORMAT! ///T40, ' APPENDI X 

IX, • COEFFICIENT' ) 

F0RMATC3I5.F5. 

FORMAT ( ////AX, 

FORMAT (////4X, 

FORMAT! 7F10.1) 

FORM AT ( F 1 0. 3 , F 1 0 . 1 , 2F 1 0 . 2 , 1  5 ) 
FORMAT (//5X, 'MISSILE    LENGTH    = 

■(L/D    =     '  ,F3.1  ,' )  '//6X,  • ADA' 

«RPM«  ,3X,  'P-HAT  ',  3X,  «Y', 5X, 

•2X,«CY(V+) ' , IX,' CYPOW' ) 
FORMAT ( F 10. 1, F 1 0. 2, 4F 10.1) 
FORMAT(«0',4X,F4.1,1X,1PE12.5,2X,OPF5.1,2X,F6.0,2X, 

F4.1,2X.F4.2,7(2X,F5.2) ) 
FORMAT! 7F10.3) 


• ,F6.3, IX, • INCHES' , IX, 
6X, ' RE(L) ' ,7X,' V ,5X, 
CY' ,3X, 'CY(-) » ,2X,'CY(+) • 


9    STOP 
END 


FUNCTION    PIF2(X,FLIST,N,GLIST) 
SUBPROGRAM    FUNCTION    2ND    ORDER    LOCKUP 


DIMENSION  FLIS 
BLIF(P,0,R,S,T 
IF     (X-FLIST(N) 

1  I=N-1 
GO    TO    5 

2  IF(X-FLIST(D) 

4  1=1 

5  K  =  l 

GO    TO    30 

6  K  =  2 


T( 100) ,GLIST(100) 

)  =  (  (  0-P)-'(S-T)/(R-0)+S) 

)     2,1,1 


4,4,6 
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9 
30 
10 
11 

12 
13 
14 
15 

16 
17 

18 


DH    8       I=1.N 

IF    (X-FLIST(I) 

CONTINUE 

I=N 

1=1-1 

BLIF1=BLIF( X 

IF     (K-l)        11 

PIF2    =     BLIF1 

RETURN 

IFC (1+2 )-N) 

IF( (I-l)-l) 


)     9,9,8 


FLIST( 
11.  12 


I)  ,FLIST(I-H)  ,GLIST(I  )  .GLISTU  +  1)  ) 


13,  13, 16 

15,14,14 

IF    (ABS     (FLIST(  I-D-X)    -ABS    (  F  L  I  ST  (  I +  2  ) -X  )  )        16,15,15 
L=I+2 
GO    TO    17 
L=I-1 

BLIF2=ELIF(X,FLIST( I) ,  FLIST ( L ) , GL 1ST (  I),GLIST(L)  ) 
PIF2=BL IF(X,FLIST(I  +  1),FLIST(U,BLIF1,BLIF2) 
RETURN 
END 
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